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We present a way to realize a n-qubit controlled phase gate with superconducting quantum 
interference devices (SQUIDs) by coupling them to a superconducting resonator. In this proposal, 
the two logical states of a qubit are represented by the two lowest levels of a SQUID. An intermediate 
level of each SQUID is utilized to facilitate coherent control and manipulation of quantum states 
of the qubits. It is interesting to note that a n-qubit controlled phase gate can be achieved with n 
SQUIDs by successively applying a n/2 Jaynes-Cummings pulse to each of the n — 1 control SQUIDs 
j^-^ | before and after a n Jaynes-Cummings pulse on the target SQUID. 
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Quantum computing has attracted much interest since it was clear that quantum computers are in principle able 
to solve hard computational problems much more efficiently than classical computers [1-3]. In the past decade, 
various physical systems have been considered for building up quantum information processors. Among them, cavity 
, QED analogs with solid state systems are particularly appealing. Theoretically, it was predicted earlier that the 
' strong coupling limit, which is difficult to achieve with atoms in a microwave cavity, can readily be realized with 
superconducting charge qubits [4,5], superconducting flux qubits [6], or semiconducting quantum dots [7]. Recently, 
the strong coupling cavity QED has been experimentally demonstrated with superconducting charge qubits and flux 
qubits [8,9] and semiconductor quantum dots embedded in a microcavity [10-12]. The results of these experiments 
make solid state qubit cavity QED a very attractive approach to quantum information process. 
, It is known that a quantum computing network can be decomposed into two-qubit gates and one-qubit rotations 
[13]. So far, a large number of theoretical proposals for realizing two-qubit gates have been presented with many 
physical systems. Moreover, two-qubit CNOT gates or controlled phase gates have been experimentally demonstrated 
pi ' in cavity QED [14], ion trap [15], NMR [16], quantum dots [17], and superconducting charge qubits [18]. On the other 
hand, research on quantum computing has recently moved towards the physical realization of multiqubit quantum 
gates. Several schemes for realizing three-qubit Toffoli gates have been proposed with neutral atoms in an optical 
lattice [19] or hybrid atom-photon qubits via cavity QED [20]. In addition, experimental realization of a controlled 
phase gate in a three qubit NMR quantum system has been reported recently [21]. 

In this paper, our goal is to present a way to realize a n-qubit controlled phase gate with superconducting quantum 
interference devices (SQUIDs), coupled to a superconducting resonator. In this proposal, the two lowest levels |0) 
and |1) of each SQUID represent the two logical states of a qubit while a higher-energy level |2) of each SQUID is 
used to facilitate coherent control and manipulation of quantum states of the qubits (Fig. 1). The method presented 
here essentially operates by having the resonator to be: (a) resonant with the |0) <-» |2) transition of each SQUID 
sequentially, while (b) largely detuned from the |1) <-> |2) transition and the |0) «-> |1) transition of each SQUID during 
the gate operation. We find that a n-qubit controlled phase gate can be achieved with n SQUIDs, by successively 
applying a 7r/2 Jaynes-Cummings pulse to each of the n — 1 control SQUIDs before and after a n Jaynes-Cummings 
pulse on the remaining SQUID (the target qubit). 

As shown below, this scheme has the following advantages: (i) No tunneling between the qubit levels |0) and |1) is 
required so that the storage time of each qubit can be made very long; (ii) No measurement on SQUIDs or photons 
is needed, therefore the operation is simplified; (iii) No auxiliary SQUIDs are needed, thus saving precious hardware 
resources; and (iv) Coupling constants of each SQUID with the resonator could be non-identical to accommodate 
inevitable nonuniformity in device parameters. 

The motivation for this work is threefold: (i) SQUIDs have recently attracted much attention in quantum infor- 
mation community. Reasons for this are that SQUIDs are relatively easy to scale up and have been demonstrated 
to have relatively long decoherence times among solid state qubits [9,22-27] and therefore have been considered as a 
promising candidate for building up superconducting quantum computers/information processors [28-37]. (ii) In the 
present proposal, the strong coupling between a SQUID and the resonator is analogous to atomic cavity quantum 
electrodynamics (QED). Although many schemes for realizing two-qubit quantum gates with a variety of physical 
systems using cavity QED [5,32-34,38-45] or trapped ions [46-49] have been proposed, how to realize a multiqubit 
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controlled quantum gate based on cavity QED or trapped ions has not been thoroughly investigated, (iii) As is 
well known, multiqubit controlled quantum gates are of great importance to constructing quantum computational 
networks, realizing quantum error correction protocols, and implementing quantum algorithms. 

This paper is organized as follows. In Sec. II, we review the basic theory of a SQUID coupled to a single-mode 
resonator or driven by a classical microwave pulse. In Sec. Ill, we show how to realize a three-qubit controlled phase 
gate with three SQUIDs in a resonator. We then discuss how the method can be generalized to implement a n-qubit 
controlled phase gate with n SQUIDs in a resonator. In Sec. IV, we give a discussion on the experimental issues. A 
concluding summary is given in Sec. V. 

II. BASIC THEORY 

The SQUIDs considered throughout this paper are rf SQUIDs each consisting of a Josephson tunnel junction 
enclosed by a superconducting loop (the size of an rf SQUID is on the order of 10 ^m — 100 /xm). The Hamiltonian 
for an rf SQUID (with junction capacitance C and loop inductance L) has the usual form [25] 

where the magnetic flux threading the ring, and Q, the total charge on the capacitor, are the conjugate variables 
of the system (with the commutation relation [<j>, Q] — ih), <fr x is the static (or quasistatic) external magnetic flux 
applied to the ring, and Ej = I c $ /2tt is the Josephson coupling energy (I c is the critical current of the junction and 
$0 = h/2e is the flux quantum). 

A. SQUID coupled to a single-mode resonator 

Consider a SQUID coupled to a single-mode resonator. The SQUID is biased properly to have the A-type three 
lowest levels, which are denoted by |0) , |1) and |2), respectively (Fig. 1). Suppose that the coupling of |0) , |1) and 
1 2) with other levels of the SQUID via the resonator is negligible, which can be readily achieved by adjusting the 
level spacings of the SQUID. It can be shown that when the resonator is resonant with the |0) <-» |2) transition while 
far-off resonant with the |1) <-> |2) transition and the |0) <-> |1) transition of the SQUID, the interaction Hamiltonian 
in the interaction picture, after making the rotating-wave approximation, can be written as [32] 

#/ = ?i(.ga + |0)(2|+h.c.). (2) 

Here, a + and a are the creation and annihilation operators of the resonator mode with frequency v r , and g is the 
coupling constant between the resonator mode and the |0) <-> |2) transition of the SQUID. For a superconducting 
one-dimensional standing-wave resonator, the expression of g is given by 

where M sr is the mutual inductance between the SQUID and the resonator, Lq is the inductance per unit length of 
the resonator, I is the length of the resonator, v r is the frequency of the resonator mode with wavelength A, and x is 
the center of the SQUID in the resonator. 

The Hamiltonian (2) actually has the same form as the Jaynes-Cummings Hamiltonian of a 2-level system resonant 
with a single- mode cavity field (a Hamiltonian well known in quantum optics). In the case when the resonator is 
initially in the photon- number state \n), the time evolution of the states of the system, governed by the Hamiltonian 
(2), is described by 

|0) \n) — > cos \fngt |0) \n) — i sin \fngt |2) \n — 1) , 

|2)|n) -> -ism\/n + lgt |0) \n + 1) + cos + lgt |2) \n) . (4) 

Note that the coupling strength g may vary with different SQUIDs due to non-uniform device parameters and/or 
non-exact placement of SQUIDs in the cavity. Therefore, hereafter we replace g by g\, g 2 , and g n for SQUIDs 
1,2,..., and n, respectively. 

B. SQUID driven by a classical microwave pulse 
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Now let us consider a SQUID driven by a classical microwave pulse with the magnetic component B /JU) (r,t) = 
B MU ,(r) cos (2Tru^ w t + 4>) . Here, B /JU) (r), u^ w , and <fi are the magnetic field amplitude, frequency, and phase of the 
microwave pulse. It can be shown that if the microwave pulse is resonant with the |1) <-» |2) transition but far-off 
resonant with the |0) <-> |2) transition and the |0) <-» |1) transition of the SQUID, then the interaction Hamiltonian 
in the interaction picture is given by 

ff/ = !(lW |l)(2|+h.c), (5) 
where Q12 is the Rabi frequency of the pulse, which takes the following form [32] 

n 12 (t) = -^(l|$|2)^B^(r)-c/S. (6) 
From the Hamiltonian (5), it is straightforward to see that a pulse of duration t results in the following rotation 

|1) ^cos^|l)-ze-^sin^|2>, 

|2) ->-i e *sin^i|l)+cos^i|2>. (7) 



III. MULTIQUBIT CONTROLLED PHASE GATE WITH SQUIDS 

In this section, for clarity, we will first give an explicit description on how to realize a three-qubit controlled phase 
gate with three SQUIDs coupled to a microwave resonator. We will then discuss how to extend the method to obtain 
a n-qubit controlled phase gate with a larger number n. 

A. Three-qubit controlled phase gate 

For three qubits, there are a total number of eight (2 3 ) computational basis states, denoted by |000) , |001) , | 111) , 
respectively. A three-qubit controlled phase gate is described by 

1 000) -» |000) , 1 100) -» 1 100) , 

|001) -» |001) , 1 101) -» |101) , 

|010) -» |010) , 1 110) -» |110) , 

|011> — |011> , |111> — - |111> , (8) 

which implies that if and only if the two control qubits (the first two qubits) are in the state |1), a phase flip happens 
to the state |1) of the target qubit (the last qubit) and nothing happens otherwise. 

To realize this gate, consider SQUIDs (1,2,3) each having the A-type level configuration as depicted in Fig. 1. 
The transition between any two levels for each SQUID is initially far-off resonant with the resonator (e.g., via prior 
adjustment of the level spacings) and the cavity mode is initially in the vacuum state |0) c . 

The operations for realizing the three-qubit controlled phase gate are listed below: 

Step (i): Apply a tt microwave pulse (Qi2T^ w — n, where 77™ is the pulse duration) with cf> = —tt/2 to SQUID 
1 [Fig. 2(a)]. The pulse is resonant with the |1) <-> |2) transition of SQUID 1. After the pulse, the transformation 
|1> — > j2> of SQUID 1 is obtained. 

Step (ii): Bring the |0) <-> |2) transition of SQUID 1 to resonance with the resonator for an interaction time 
n = tt/ (2 5l ) [Fig. 2(b)], resulting in |2> x |0) c - -i \0) 1 |l) c . 

Step (iii): Bring the |0) «-> |2) transition of SQUID 2 to resonance with the resonator for an interaction time 
T 2 = tt/ (252) [Fig. 2(c)]. As a result, the states |0) 2 |0) c , |1) 2 |0) c , and |1) 2 |l) c remain unchanged, while the state 
|0) 2 |l) c changes to -i |2) 2 |0) c . 

Step (iv): Bring the |0) <-» |2) transition of SQUID 3 to resonance with the resonator for an interaction time 
t 3 = 7r/g 3 [Fig. 2(d)], resulting in |0) 3 |l) c — > — |0) 3 |l) c while no change for the states |0) 3 |0) c , |l) 3 |0) c , and 

|l) 3 |l)c- 

Step (v): Bring the |0) <-> |2) transition of SQUID 2 to resonance with the resonator for an interaction time 
T 2 = tt/ (252) [Fig. 2(c)]. As a result, the states |0) 2 |0) c , |1) 2 |0) c , and |1) 2 |l) c remain unchanged, while the state 
|2) 2 |0) c becomes -i|0) 2 |l) c . 
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Step (vi): Bring the |0) <-> |2) transition of SQUID 1 to resonance with the resonator for an interaction time 
n = tt/ (2 5l ) [Fig. 2(b)], resulting in [0^ |l) c - -i \2) 1 |0) c . 

Step (vii): Apply a n microwave pulse with = tt/2 to SQUID 1 [Fig. 2(a)]. The pulse is resonant with the 
|1) <-> |2) transition of SQUID 1. After the pulse, the transformation |2) — ► |1) of SQUID 1 is achieved. 

The states of the whole system after each step of the above operations are summarized below: 

|100)|0) c |200)|0) c -i|000)|l) c -|020)|0) c 

|101)|0) c stop(i) |201)|0) c stop(ii) -i|001)|l) c stcp(iii) -|021)|0) c 

|110) |0) c ~ * |210)|0) c ~* -t|010)|l) c ~* -*|010)|l) c 

|111) |0) c |211) |0> c -i|011)|l) c -<|011>|l) c 



-|020)|0) c i|000)|l) c |200) |0) c |100) |0) c 

stcp(iv) -|021)|0) c stcp(v) *|001)|l) c stcp(vi) |201)|0) c stcp(vii) |101)|0) c . . 

^ i|010)|l) c ^ i|010)|l) c ^ |210) |0) c ^ |110) |0) c ' W 

-i|011)|l) c -i|011)|l) c -|211)|0) c - jlll> |0> c 

where \ijk) is abbreviation of the state \i} 1 \ j} 2 \k) 3 of SQUIDs (1, 2, 3) with i, j, k e {0, 1, 2}. 
On the other hand, it is obvious that the following states of the system 

|000) |0) c , |001) |0) c , |010) |0) c , |011) |0) c (10) 

remain unchanged during the operation. This is because: (a) the state |0) of SQUID 1 was not effected by the applied 
microwave pulse, since the |0) <-> |2) transition and the |0) <-» |1) transition of SQUID 1 are far-off resonant with the 
applied microwave pulse; and (b) no photon was emitted to the resonator during the Step (ii) operation, when SQUID 
1 is initially in the state |0). Hence, it can be concluded from Eq. (9) that the three-qubit controlled phase gate (8) 
was achieved with three SQUIDs (i.e., the control SQUIDs 1 and 2, as well as the target SQUID 3) after the above 
process. 



B. n-qubit controlled phase gate 

A quantum controlled phase gate of n qubits (1, 2, n) is defined by the following transformation 

\hi2-in) -» (-i) llXt2 --- X4 " | i l i 2 ...ii...i n ) , (11) 

where the subscript I represents qubit I, i\ e {0, 1}, and \i\i2-..ii...i n ) is a n-qubit computational basis state. For n 
qubits, there are a total number of 2" computational basis states, which form a set of complete orthogonal bases in 
a 2"-dimensional Hilbert space of the n qubits. Eq. (11) shows that only when the n — 1 control qubits (the first 
n — 1 qubits) are all in the state |1), the state |1) of the target qubit (the last qubit) undergoes a phase flip, i.e., 

1 1 1 1) — ► — 1 1 1 1) . While nothing happens to all other 2" — 1 computational basis states. In the following, we will 

discuss how this gate can be achieved with n SQUIDs coupled to a resonator. 

The n SQUIDs are labeled by 1,2,..., and n. The first n — 1 SQUIDs (l,2,...,n— 1) act as control qubits while 
the SQUID n is the target qubit. Suppose that the n SQUIDs (1,2, ...,n) are initially decoupled from the resonator 
(initially in the vacuum state). Examining the above operations for the three-qubit controlled phase gate carefully, 
we find that the n-qubit controlled phase gate (11) can be obtained with the resonator mode returning to the original 
vacuum state, by the sequence of operators 



u = u+ ® TT u lr I ® u% r <g> I TT Uir ] ® Ui, (12) 




where n»=i Uir = Uk r U^^iy r ■ ■ ■ U2 r U\ r \ Ui and U± denote the operators on SQUID 1 with the following matrixes of 

(13) 



UtUi =1, Ux = ( _°! I ) 
= (1, 

(I = 1, 2, ...,n — 1), represented by the matrix 



in the basis states |1) 1 = (0, 1) T and \2) 1 = (1,0) T ; Ui r is the joint operator on the SQUID I and the resonator mode 
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in the basis states |0) ; |1) = (0, 1) T and |2). |0) = (1,0) T ; and U nr is the joint operator on the SQUID n and the 
resonator mode with the same matrix as the one described by (14) [in the basis states \0) n |l) c = (0, 1) T and \2) n |0) c = 

(i,o) T ]. 

From the description in the previous subsection, it can be seen that: 

(i) U\ (Ui) denotes the application of a tt microwave pulse (f2i2T Mlt) — n, where t^ w is the pulse duration) with 
4> = —it/2 (it/2) and v^ w — v\i (the |1) <-> \2) transition frequency of SQUID 1) to SQUID 1; 

(ii) Uir corresponds to the operation of bringing the |0) <-> |2) transition of SQUID I (I = l,2...,n— 1) to resonance 
with the resonator for an interaction time ti = tt / (2gi) (i.e., a it /2 Jaynes-Commings pulse); and 

(hi) U% r indicates two it/2 Jaynes-Commings pulses resonant with the |0) <-> |2) transition of SQUID n, which are, 
when combined together, equivalent to a tt Jaynes-Commings pulse (i.e., r„ = tt / g n ). 

We remark that all other basis states of the system involved in each step of the above operations, which form a 
complete set of orthonormal states together with the basis states described above, are not affected by (a) setting the 
microwave pulse to be largely detuned from the |0) <-» |2) transition and the |0) <-» |1) transition of SQUID 1 and (b) 
setting the |1) <-> |2) transition and the |0) <-> |1) transition of each SQUID far-off resonant from the resonator mode. 

The method presented here for realizing the n-qubit controlled phase gate (11) is the extension of the three-qubit 
version (8) described in the previous subsection. This can be seen as follows. From Eq. (9), one can see that the 
three-qubit controlled phase gate (8) was realized essentially through the operation of step (iv). This operation lead 
to a phase shift on the state of SQUID 3 (i.e., |0) — > — |0) and |1) — ► |1)) with the aid of the photon, when the two 
control SQUIDs 1 and 2 are initially in the basis state |11) . But, when the two control SQUIDs 1 and 2 are initially 
in the other basis states, this operation results in no change to the state of SQUID 3, due to the fact that no photon 
was left in the resonator mode after the step (iii). Similarly, the realization of the n-qubit controlled phase gate (11) 
described above is mainly based on the operation described by U% r , which causes a phase shift on the state of SQUID 
n (the target) with the assistance of the photon when the n — 1 control SQUIDs (1, 2, ...,n — 1) are initially in the 
basis state |ll...l) . Note that only for this initial basis state, the photon originally created by the operation Ui r Ui 
would remain in the resonator mode after the operation Y\i=2 Uir- I n addition, similar to the operations of steps 
(v)-(vii) for the three-qubit controlled phase gate (8), the operation described by {/+ (g> ^nl=n-i Uirj has the photon 
emitted back into the resonator mode from the system of SQUIDs (2, 3, ...n — 1) and finally absorbed by SQUID 1. 
Hence, the resonator mode returns to the original vacuum state and the n SQUIDs (1, 2, n) are back to the initial 
states except a phase flip to the n-qubit basis state |ll...l). 

Before closing this section, some points may need to be addressed here. 

(a) . The irrelevant SQUIDs in each step of the operation need to be decoupled from the resonator/pulse during 
the resonator /pulse- SQUID interaction. The resonator mode needs to be not excited during the application of the 
microwave pulse. In addition, the coupling of the levels |0), |1) , and |2) with the other levels should be negligible for 
each SQUID. In principle, these conditions can be satisfied by adjusting the level spacings of the SQUIDs. Note that 
for a SQUID, the level spacings can be changed readily by varying the external flux <fr x or the critical current I c (e.g., 
for variable barrier rf SQUIDs) [50]. 

(b) . It is not necessary to have a single-mode resonator since for a multi-mode resonator one can choose one mode 
to interact with the SQUIDs while have all other modes well decoupled from the three lowest levels of the SQUIDs 
(e.g., with proper device parameters). 

(c) . The method presented here is applicable to a ID, 2D, or 3D microwave resonator /cavity as long as the conditions 
described above can be met. 

(d) . As is well known, a n-qubit controlled-NOT (CNOT) gate (known as the Toffoli gate for n = 3) can be obtained 
by combining the n-qubit controlled phase gate (11) with two single-qubit Hadamard gates, which are performed on 
the target qubit before and after the n-qubit controlled phase gate (11) respectively (see Fig. 3). Each of the single- 
qubit Hadamard transformations, |0) — + ^- (|0) + |1)) and |1) — > (|0) — |1)) , can be done with a tt/2 microwave 
pulse resonant with the |0) <-> |1) transition of the target SQUID qubit. When combining with the above quantum 
controlled phase gate operations, one obtains the n-SQUID qubit CNOT gate. 

(e) . The present method provides a simple way to implement a multi-qubit CNOT gate with SQUIDs coupled to 
a resonator. To see this, let us consider the simple case of n = 3. It is known that construction of a Toffoli gate 
requires at least six two-qubit CNOT gates and ten single-qubit gates (i.e., two Hadamard, one phase, and seven 7r/8 
gates) [51]. Note that a two-qubit CNOT gate consists of a two-qubit controlled phase gate and two single-qubit 
Hadamard gates as described above. Therefore, using the conventional gate-constructing technique, at least 28 steps 
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of operations will be necessary to realize a Toffoli gate, assuming that the realization of a single-qubit gate or a 
two-qubit controlled phase gate requires one-step operation only. However, as discussed above, the present method 
only needs 9 steps of operations. That is, seven steps of operations for the three-qubit controlled phase gate (8) plus 
two steps of operations for the two Hadamard gates. Finally, it is obvious that the simplicity of the present method 
in constructing a n-qubit CNOT gate may become more apparent with the increment of n, when compared with the 
use of the conventional gate-decomposition protocol. 



In this section we discuss issues that are important to experimental implementation. Without loss of generality, let 
us consider performing a n-qubit controlled phase gate with n identical SQUIDs (1, 2, n) at locations where the B r 
fields are the same (e.g., antinodes of the cavity field). In this case, we have gi = g (I = 1, 2, n). For the method 
to work, the total operation time r = 2n (r r + r a ) + 2r fiw [r r = ir/ (2g)] should be much shorter than the energy 
relaxation time j^ 1 of the level |2) , and the lifetime of the resonator mode = Q/2ttv t , where Q is the (loaded) 
quality factor of the resonator. Here, r a is the typical time required for adjusting the level spacings of a single SQUID. 
In addition, direct coupling between SQUIDs needs to be negligible since this interaction is not intended. 

In principle, these requirements can be realized, since: (i) r r can be reduced by increasing the coupling constant 
g, (ii) r a can be shortened by rapid adjustment of the level spacings of the SQUIDs, (iii) can be increased by 
employing a high-Q resonator so that the resonator dissipation is negligible during the operation, (iv) the SQUIDs 
can be designed so that the energy relaxation time j^ 1 of the level |2) is sufficiently long, and (v) direct interaction 
between SQUIDs is negligible as long as the following condition can be met 



where H s ^ s is the interaction energy between the two nearest SQUIDs, H s - r is the SQUID-resonator interaction 
energy, and the SQUID-microwave interaction energy. 

It is straightforward to see that the condition (15) can be realized if 



Here, M is the mutual inductance between the two nearest SQUIDs I and I + 1 (I = 1,2, ...,n — 1), <pf) {4^^) 
= (i\ $ |j) /$o is the magnetic dipolc coupling matrix element between levels |z) and \j) of SQUID I (I + 1), and 
ij = 01,02, or 12. 

For the sake of definitiveness, let us consider the experimental feasibility of realizing a three-qubit controlled phase 
gate using SQUIDs with the parameters listed in Table 1. Note that SQUIDs with these parameters are readily 
available at the present time [22,23,52]. With the choice of these parameters, the SQUIDs have the desired three-level 
structure as depicted in Fig. 1. For a superconducting one dimensional standing- wave CPW (coplanar waveguide) 
resonator with the parameters listed in Table 1 and SQUIDs placed along the cavity axis (Fig. 4), one has M sr <~ 100 
pH. When each SQUID is located at one of the antinotes of the resonator mode (Fig. 4), a simple calculation shows 
g <~ 7.5 x 10 9 s _1 , resulting in r r w 0.2 ns. With the choice of t^ w ~ r a <~ r r , one has r ~ 2.8 ns, which is much smaller 
than ~ 3.2 ^s and kT x <~ 41.7 ns for a resonator with Q r = 3 x 10 3 . Note that superconducting CPW resonators 
with higher quality factors have been demonstrated by recent cavity QED experiments with superconducting charge 
qubits [8]. 

For a resonator with v r = 11.4 GHz, the wavelength of the resonator mode is A ~ 10.5 mm. When each SQUID is 
placed at an antinode of the B r field (Fig. 4), one has D ~ 5.3 mm, where D is the distance between the two nearest 
SQUIDs. A simple estimate gives M < 0.1 aH, resulting in ( < 10~ 8 for the parameters considered above. Thus, the 
condition of negligible direct coupling between SQUIDs is very well satisfied. 

The above analysis demonstrates that the realization of a three-qubit controlled phase gate is possible using SQUIDs 
and a resonator. We remark that a quantum controlled phase gate with a larger number of qubits can in principle 
be obtained by increasing the length of the resonator though the conditions of r <C j^ 1 , kT x becomes increasingly 
difficult to satisfy. 

We emphasize that the primary purpose of this work is to provide a new approach to implement a multiqubit 
quantum controlled phase gate with SQUIDs. However, we note that: (a) when coupled to a cavity mode, many 
physical qubit systems (such as atoms, quantum dots, and superconducting charge qubits) have the same type of 
qubit-cavity interaction described by the Hamiltonian (2), and (b) the condition, i.e., the |0) <-> |2) transition being 
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resonant while the transition between any other two levels being far-off resonant with the resonator, can always 
be obtained via the adjustment of the level spacings (e.g., for quantum dots and atoms, the level spacings can be 
changed via adjusting the external electric field [53]). Therefore, it is straightforward to show that the method can 
be generalized to realize the multiqubit controlled phase gate in other types of qubit systems with the A-type three 
level configuration within cavity QED. 

V. CONCLUSION 

Before conclusion, we should point out that the idea of tuning the individual qubits in and out of resonant from 
the cavity mode was previously proposed for superconducting charge qubits [54]. Rather, our scheme is for a different 
system and it differs in the details of both the qubits and the coupling structure. In our case, we consider a system 
consisting of flux qubits (SQUIDs) coupled to a microwave resonator/cavity, while the system described in [54] 
comprises charge qubits and a LC-oscillator mode in the circuit. In addition, the idea of performing a phase shift 
via the assistance of the cavity photon was proposed earlier for realizing a two- qubit quantum controlled phase gate 
with trapped ions [46]. However, to the best of our knowledge, our scheme is the first to demonstrate that a quantum 
controlled phase gate with a large number of qubits can in principle be achieved within cavity QED which is of great 
importance. 

In summary, we have presented a method to realize a multiqubit controlled phase gate with SQUIDs coupled to 
a microwave resonator, which operates essentially by exchanging a single photon between the control SQUIDs and 
the resonator mode before and after a phase shift performed on the target SQUID. The method has the following 
advantages, (i) Only one SQUID interacts with the microwave pulse; (ii) No auxiliary SQUIDs or measurement is 
needed during the entire operation, thus the hardware resources is significantly reduced and the operation is greatly 
simplified; (iii) As tunneling between the qubit levels |0) and |1) is not required during the operation, decay from the 
level |1) can be made negligibly small during the operation (via prior adjustment of the potential barrier between 
the qubit levels |0) and |1) [50]) and therefore the storage time of each qubit can be made much longer; (iv) The 
coupling constants of SQUIDs with the resonator could be different, which makes the present proposal much easier 
to implement since neither identical SQUIDs nor exact placement of SQUIDs is needed; and (v) The method can in 
principle be applied to obtain a n-qubit controlled phase gate with a large number n. Finally, as discussed above, 
the present method is quite general, which can be applied to implement a multiqubit controlled phase gate for other 
types of physical qubit systems with the A-type three- level structure within cavity QED. 
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Figure and Table Captions 



FIG. 1. Level diagram of a SQUID with the three lowest levels |0) , |1) and |2) forming a A-type structure. 

FIG. 2. Illustration for the change of the level structure (reduced) of three SQUIDs (1,2,3) during a three-qubit 
controlled phase gate performance. In (a), (b), (c), and (d), figures from left to right represent the level structures 
for SQUIDs 1,2, and 3, respectively; the non-identical level spacings of the SQUIDs could be caused by nonuniform 
device parameters. g\, g 2 , and gs are the resonant coupling constants between the resonator mode and the |0) <-» |2) 
transition of SQUIDs 1,2, and 3, respectively. The difference for gi, g 2 , and g% is due to device parameter non- 
uniformity or non-exact placement of each SQUID. u^ w is the frequency of the applied microwave pulse while v\2 is 
the |1) <-> |2) transition frequency for SQUID 1. In (a), the level spacings of SQUID 1 is set to be much different from 
that of SQUIDs 2 and 3, such that SQUIDs 2 and 3 are decoupled from the applied pulse. The transition between 
any two levels linked by a dashed line is far-off resonant with the resonator mode. 

FIG. 3. Relationship between a n-qubit controlled NOT gate and a n-qubit controlled phase gate. For the circuit 
on the left side, the element denoted by corresponds to a controlled NOT (with n — 1 controls on the filled circles); 
if the n — 1 controls are all in the state |1) , then the state at © is bit-flipped. On the other hand, for the circuit 
on the right side, the element Z represents a Pauli rotation a z , i.e., a phase flip operation (with n — 1 controls on 
the filled circles). Namely, if the n — 1 control qubits are all in the state |1) , then the state |1) at Z is phase- flipped 
as |1) — > — |1) while nothing happens to the state |0) at Z. In addition, the element containing H corresponds to a 
Hadamard transformation described by |0) — > -j= (|0) + |1)) while |1) — > (|0) — |1)) . 

FIG. 4. Sketch of the setup for three SQUIDs (1,2,3) and a standing-wave quasi-one dimensional CPW resonator 
(Not drawn to scale). Each SQUID is placed in the plane of the resonator between the two lateral ground planes (i.e., 
the x-y plane) and at an antinode of the B r field. The two curved lines represent the standing-wave B r field, which 
is in the ^-direction. 

TABLE 1. Parameters for a SQUID-resonator. (3l is the SQUID's potential shape parameter, R is the SQUID's 
effective damping resistance, and S is the surface bounded by the loop of the SQUID with width a and length b. r ) 2 1 
(7f ) is the energy relaxation time of the level |2) (|1)). ^02 (^12) is the |0) <-> |2) (|1) <-> |2)) transition frequency. 
4>ij = (i\ $ \ j) /$ is the magnetic dipole coupling matrix clement between levels \i) and \j) (i = 0, 1; j = 2). I is the 
length of the quasi-one dimensional CPW resonator, A is the wavelength of the resonator mode with frequency v r , d 
is the gap between the center conductor and the adjacent ground plane, w is the width of the center conductor, t is 
the width of each ground plane, L is the inductance per unit length of the waveguide, and e e is the effective relative 
dielectric constant. 
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